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Using polarized neutron scattering we establish that the magnetic order in Lai.4sNdo.4Sro.i2Cu04 
is either (i) one dimensionally modulated and collinear, consistent with the stripe model or (ii) two 
dimensionally modulated with a novel noncollinear structure. The measurements rule out a number 
of alternative models characterized by 2D electronic order or ID helical spin order. The low-energy 
spin excitations are found to be primarily transversely polarized relative to the stripe ordered state, 
consistent with conventional spin waves. 

PACS numbers: 74.72.Dn, 75.30.Fv, 75.50.Ee, 75.70.Kw 



One of the most striking and robust features in the 
phenomenology of hole-doped copper oxide supercon- 
ductors is the four-fold incommensurate (IC) pattern of 
magnetic neutron scattering peaks centered on the an- 
tiferromagnetic (AFM ) wave vector of the square CuC>2 
lattice. This pattern is found in the magnetic excita- 
tion spectrum of YBa2Cu306+j, and La2- a; (Sr,Ba) 2 ;Cu04 
over wide doping ranges jl]. Near x — 1/8 certain La- 
based materials develop an elastic IC magnetic compo- 
nent accompanied by second order harmonics around the 
structural Bragg peaks @, 0, 0] . One school of thought as- 
sociates these features with one-dimensional (ID) charge 
modulations separating AFM antiphase bands on the 
Cu02 layers In this 'stripe' model the four- fold pat- 
tern is a superposition of two two-fold patterns, aris- 
ing from spatially separated stripe domains, each with 
charge modulations along one of the two Cu-0 bond di- 
rections. Static stripes, posited to occur near x = 1/8, 
are thought to compete with superconductivity [6(, but 
dynamic stripes could play a role in the formation of the 
superconducting state [7|. 

Recently, however, the stripe picture has been called 
into question. Several new experimental findings point 
to the existence of 2D charge density wave order in the 
ground state of hole-doped cuprates [1, Q . In addition, 
the dimensionality of the spin excitation spectrum is a 
subject of debate [13]. Furthermore, an increasing num- 
ber of phases exhibiting novel 2D electronic order have 
been explored theoretically [ill ] , including orbital current 
correlations [Hj], checkerboard-type orderings of Cooper 
pairs |l3j, and 2D diagonal stripes [14j . The possibility 
that features previously attributed to stripes might be 
signatures of a more elaborate 2D ordering makes it vital 
to obtain information with techniques that can separate 



different spatial arrangements of spin and charge. 

Here we report a study of the magnetic order and dy- 
namics in La 1 . 6 _ a: Ndo.4Sr x Cu04 (LNSCO; x = 0.12) by 
polarized-neutron scattering. Substitution of Nd for La 
stabilizes a low-temperature tetragonal structure, which 
permits the formation, below ~ 50 K, of a robust spin- 
charge-ordered phase with suppressed superconducting 
transition temperature The relatively large ordered 
Cu moment ~ 0.10/zs [3| allows detailed neutron po- 
larisation analysis. Below Tn<i — 3K, Nd-Cu coupling 
causes alignment of the Nd spins along the c axis with 
the same ordering vectors as the Cu spins Q. At tem- 
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FIG. 1: Models for the magnetic order in La-based cuprates 
at x = 1/8 doping, (a) One-q domain with charge stripes 
(lines of open circles) and collinear spin order [sj. (b) One-q 
domain with helical spin order, (c) Collinear two-q structure 
of a diagonally modulated commensurate AFM [lj]. (d) Two- 
q order of charge and spins in a noncollinear structure. 
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FIG. 2: (a) Reciprocal space of the square Cu02 lattice show- 
ing the three quartets of magnetic peaks investigated. The 
peaks are displaced from AFM wave vectors by ±qi = (±5, 0) 
and ±q2 = (0, ±<5) with 8 ~ x = 0.12. Also shown are the 
axes to which the polarisation P and magnetic scattering are 
referred, (b)-(c) Unpolarized elastic scans through the ex- 
pected Bragg peak positions for a 7r-spiral [ljj. 



peratures well above TNd the Nd order is unlikely to in- 
fluence the ordering pattern and in-plane direction of the 
Cu spins. We therefore believe that our study of static 
and dynamic properties of LNSCO at 10 K has direct 
relevance to the magnetic behaviour of Nd-frec higher 
temperature superconductors. Our results are most nat- 
urally understood in terms of a ID modulation of the 
AFM order, consistent with the occurrence of stripes, al- 
though an exotic 2D noncollinear order is also possible. 

Fig. [T] shows four models yielding magnetic diffrac- 
tion patterns with principal Fourier components Q = 
(1/2, 1/2) ± qi and (1/2, 1/2) ± q 2 and equivalent wave 
vectors, as found experimentally (Fig. [2]). Figure [TJ a) 
represents the conventional view [2| that the quartet is 
due to incoherent superposition of scattering from two 
equally populated domains with collinear one-q spin or- 
der and orthogonal propagation vectors. From the peak 
positions alone, this model is indistinguishable from a 
model of two domains, each with helical one-q order 15 1 
as sketched in Fig. [lib). The correct four-fold diffraction 
pattern is also produced by the collinear two-q "diagonal 
stripe" picture [3] i n Fig- He) and by the noncollinear 
two-q checkerboard structure shown in Fig.[TJd). 

Our experiment was performed on the IN20 triple- 
axis spectrometer at the ILL operated in Heusler-Heusler 
configuration with a pyrolytic graphite filter to suppress 
higher order contamination of the scattered beam. The 
crystal (T c = 6.8 K), grown by the floating-zone method, 
contained two grains separated by ~ 1°. The sample 
was mounted with the [001] axis vertical in a He cryo- 
stat. A final neutron energy of 34.8 meV allowed ac- 
cess to IC quartets surrounding several equivalent AFM 
wave vectors in the (H, K, 0) reciprocal lattice plane, de- 
noted (H, K) for short. The polarisation vector P of the 
neutron beam at the sample position was oriented along 



the x, y and z directions (see below and Fig. [2]). The 
scattered neutrons were recorded in spin-flip (SF) and 
non-spin- flip (NSF) channels, according to whether their 
spins had flipped or not on scattering. Corrections for 
the measured beam polarisation P = 0.86 were applied. 
Most measurements were performed at 10 K > T c > TNd, 
but some data taken at 1.7 K were used to fix the Bragg 
peak line shapes from the Nd magnetic order. 

The cross section for scattering of polarized neutrons 
consists of a purely nuclear term, a purely magnetic term 
and a nuclear-magnetic interference term [la ] . We as- 
sume that the latter can be neglected at the IC wave 
vectors of interest, and that the nuclear spins are unpo- 
larized. The magnetic term contains two features that 
make it possible to determine electronic spin directions. 
First, magnetic scattering originates only from electronic 
spin components Sj_ perpendicular to Q. Second, SF 
scattering is caused by spin correlations perpendicular to 
P [3]. We denote coherent nuclear scattering by N and 
magnetic scattering by M x , M y and M z , where M a is 
proportional to the time Fourier transform of the corre- 
lation function {S a (— Q, 0)5 a (Q, t)) [I^]. We define axes 
such that x is parallel to Q, and y and z are perpendicu- 
lar to Q in and out of the scattering plane, respectively, 
(Fig. [2|). With these axes, M x is identically zero, and in 
the absence of a single-domain chiral structure [18| the 
intensities in the SF and NSF channels with P parallel 

M y + M z + B SF , 

NSF = N + ^NSF, 

-M z +i?NSF, where 
-Bsf and -Bnsf are the backgrounds. These expressions 
apply both to elastic (t = oo) and inelastic scattering. 

We first treat elastic scattering, in which case M y , M z 
are proportional to the squares of the corresponding or- 
dered spin components. At each of the three IC quartets 
indicated in Fig. [2] the neutron count rate was recorded 
in some or all six polarisation channels at each tempera- 
ture. The purpose of probing several zones is to vary the 
orientation of Q relative to the ordered spin direction so 
that the magnetic cross sections M y and M z change. 

Fig. [3] shows data obtained by scanning the sample ro- 
tation angle 9 through the (1/2, 1/2) + qx satellite peak 
with P || z. The SF and NSF channels then sense M y and 
M z , respectively. The strong NSF signal at 1.7 K [Fig. 
121a)] implies significant elastic scattering from magnetic 
moments oriented perpendicular to the Cu0 2 planes. On 
heating to 10 K the NSF signal almost vanishes [Fig. 
[3jb)] . A weak signal is present in the SF channel at both 
1.7 K and 10 K [Figs. [3jc) and (d)]. These observations 
directly confirm an earlier finding that the Nd moments 
order along the c axis @. Analysis of the data in Figs. 
[3jb) and (d) indicates that at 10 K the spin direction is 
mainly confined to the Cu0 2 planes. The data do not 
rule out a small ordered component along c, but as the 
presence or absence of such a component does not al- 
ter our conclusions, which refer to the in-plane order, we 
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FIG. 3: SF and NSF elastic scattering at (1/2, 1/2) + qi with 
P || z. The two-peak line shape is due to the presence of two 
crystallites separated by ~ 1°. The solid line in (a) is a fit to 
two Gaussians. The lines in (b)-(d) were obtained from a fit 
to the same line shape as in (a) but with the overall intensity 
scale and background allowed to vary. 



assume that at 10 K the Cu spins lie in the Cu0 2 planes. 

We now turn to the spatial arrangement and in-plane 
orientation of the Cu moments. Figure Q] shows M y at 
three IC peaks of type qi and at three peaks of type q2. 
The qi and q 2 peaks close to (1/2,1/2) have roughly 
equal intensity. By contrast, near (—1/2,3/2) the qi 
peak is clearly weaker than the q 2 peak, while the situa- 
tion is reversed near (3/2, —1/2). 

The relative intensities of the peaks in any given quar- 
tet depend on (i) the orientation of Q relative to the spin 
components contributing to the peaks and (ii) the popu- 
lation of any equivalent magnetic domains. For the one- 
q order in Fig.QJb), the qi satellites around (—1/2,3/2) 
should have same intensity as those around (3/2, —1/2), 
in disagreement with the data. This is true also for more 
than one chiral domain [l6j |. because all in-plane spin 
directions contribute equally to each peak. Considering 
next the collinear, two-q model in Fig. [He), the angle 
between Q and the unique spin direction of a single do- 
main changes only slightly between qi and q 2 in any 
given quartet, so the large intensity differences observed 
in M y near (-1/2,3/2) and (3/2,-1/2) cannot be re- 
produced by this model. Hence, we can rule out both 
in-plane helical order [Fig. [TJb)] and "diagonal stripes" 
[Fig. [He)], as well as the so-called 7T-spiral model [ijj 
because that produces q^ satellites about (1/2, l)-type 
positions which are not observed, as shown in Fig. [2Jb) 
and Etc). 

By contrast, the data are consistent with two, equally 
populated, one-q domains, each with collinear spin or- 
der, and spins (S)j of domain j approximately perpen- 
dicular to q^ as shown for one domain in Fig. [Ha). In 
this scenario, the qi and q 2 peaks near (1/2,1/2) have 
nearly identical intensities because the angle between Q 
and (S)j is the same for both domains. Conversely, near 



(—1/2,3/2), the qi (q 2 ) peak is relatively weak (strong) 
since the corresponding spins are close to being parallel 
(perpendicular) to Q. Around (3/2,-1/2) these angu- 
lar factors switch, and the qi peak should be the most 
intense, as observed. The solid lines in Fig. 0] are calcula- 
tions based on this model. As in Fig. [31 the line shape at 
each Q was fixed by fitting the large Nd signal in I§g F at 
1.7K. The corresponding curves were then multiplied by 
the modulation expected for in-plane spins (8)j perpen- 
dicular to q-,- and corrected for the form factor difference 
between the Nd-dominated 1.7 K signal and the 10 K Cu 
signal. Although the calculations do not agree in every 
detail, they clearly reproduce the salient features of the 
data. Quantitatively, a least-squares fit to all the data in 
Fig. 2] results in (S)i _L qi and (S) 2 _L q 2 with a ±3° ac- 
curacy. Our diffraction data are thus consistent with an 
incoherent superposition of orthogonal stripe domains. 

Our data are also consistent with the two-q structure 
shown in Fig. [ljd), which is a coherent superposition of 
two orthogonal stripe domains and produces the correct 
observed positions for both charge and magnetic peaks. 
This has two implications. First, one cannot infer the 
existence of ID stripes from existing measurements of 
charge and magnetic peak positions @, [f| , and second, 
although our results do not rigorously rule out a checker- 
board configuration, they do impose the noncollinear spin 
arrangement shown in Fig. [TJd) on any such model. 
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FIG. 4: Processed polarized-neutron data taken at T = 10 K, 
showing the in-plane component M y of the magnetic cross 
section at each of six IC peaks of type qi (left) and q2 (right). 
My was calculated both from the SF data (M y = 7| F — J| F ) 
and from the NSF data (M v = I^ SF — -Tnsf)- After checking 
for consistency, data thus obtained were combined. The solid 
lines represent the model calculation described in the text. 
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Wavevector Q 



Cross-section component Intensity 



(-1/2,3/2) + qi M y 7.4 ±1.3 

(-1/2,3/2) +qi M z 9.9 ±1.4 

(3/2, -1/2) - qi M y 2.4 ±1.1 

(3/2, -1/2) -qi M z 7.6 ±1.1 

Excitation components relative to [010] 

Transverse, out-of-plane 8.8 ± 0.9 

Transverse, in-plane 7.8 ± 1.4 

Longitudinal, in-plane 1.7 ±1.3 



TABLE I: Processed inelastic (Two = 5 meV) count rates per 



10 minutes at T ■ 

IsF ~ -^SF -^SF 



10 K. My and M z were obtained Irom 
7| F , respectively. Total time: 44 h. 



For inelastic scattering, M z is sensitive to spin fluc- 
tuations out of the scattering plane and M y to in-plane 
fluctuations perpendicular to Q. Table U shows the count 
rates at (-1/2, 3/2) +q x and (3/2, -1/2) - q x at hu> = 5 
meV. There is a statistically significant difference be- 
tween the values of M y at the two wave vectors. This 
shows that the low-energy spin fluctuations in LNSCO 
have a preferred direction, and it implies that they are 
not of the singlet-triplet type (for which we would expect 
My = M z ). The data are naturally explained by the 
qi stripe domain shown in Fig. [U[a). Converting from 
My to components transverse to and along (S)i shows 
(Table [TJ) that the fluctuations are predominantly trans- 
verse, consistent with spin waves. This finding supports 
theories of the cuprate spin excitation spectrum based 
on a ground state with slowly fluctuating stripelike cor- 
relations, in which the low-energy excitations resemble 
Goldstone modes of weakly coupled spin ladders [IB] ■ 

In summary, we have shown that the magnetic order in 
the spin-charge ordered cuprate Lai.4sNdo.4Sro.i2Cu04 
is either modulated in ID only (with spins perpendicular 
to the modulation direction) or takes the form of a previ- 
ously unconsidered noncollinear two-q structure. In the 
former case, it is reasonable to conclude that charge order 
is also one-dimensional, consistent with a stripe model. 
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